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Abstract 

We study, for the first time, the Schwinger mechanism for the pair production of charged scalars in 
the presence of an arbitrary time-dependent background electric field E{t) by by directly evaluating 
the path integral. We obtain an exact non-perturbative result for the probability of charged scalar 
particle-antiparticle pair production per unit time per unit volume per unit transverse momentum 
(of the particle or antiparticle) from the arbitrary time dependent electric field E{t). We find 
that the exact non-perturbative result is independent of all the time derivatives ' , where 
n = 1,2, ....oo. This result has the same functional dependence on E as the constant electric field 
E result with the replacement: E — > E{t). 
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In a classic paper in 1951 Schwinger derived the following one-loop non-perturbative 
formula 

dW e^E"^ ^ 1 

= > — e 1"=^! (1) 

d^x 47r3 ^ ^ 

for the probability of e^e~ pair production per unit time per unit volume from a constant 
electric field E via vacuum polarization [l| by using proper time method. In case of charged 
scalar field theory the corresponding result is given by 

dW _ ^ (-1)"+^ 

The pt distribution of the e"*" (or e~) production, -^ji^^, can not be obtained by using 
proper time method and a WKB approximate method was used for this purpose in Q]. A 
path integral method can also be used to obtain the px distribution [3|]. In the case of 
fermions in QED one finds B 

'Log[l-e--ni^]. (3) 



d^xd'^PT Air^ 

The corresponding result for the charged scalar production is given by 

dW \eE\ 



„2 



Log[l + e"^^^]. (4) 



d^xd^Px SvT'^ 

Recently this approach has been extended to QCD in SU(3) gauge group by directly 
evaluating the path integral and the pt distribution of the quark and gluon production 
rate are found to depend on two independent Casimir invariants: Ci = [E"'E"'] and Co = 
[dabcE°'E^E'^]^ where E"" is the constant chromo-electric field with color index a=l,2,..8 3|. 

In this paper we study, for the first time, Schwinger mechanism in the presence of an 
arbitrary time dependent electric field E{t) by directly evaluating the path integral. This is 
possible by using recently derived shift theorems for path integrals ^\ and the fact that the 
eigenvectors of the Hamiltonian for a harmonic oscillator with time dependent frequency u!{t) 
are the same as the eigenvectors of the constant frequency u problem, with the replacement 

Our result which we will derive below is the following non-perturbative formula for the 
probability of charged scalar particle-antiparticle pair production per unit time per unit 
volume per unit transverse momentum (of the particle or antiparticle) production from 
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arbitrary time dependent electric field E{t) via Schwinger mechanism 



- ' ^ 'Log[l + e ''T^W]. (5) 



This result has the remarkable feature that no time derivatives ^-^^ (where n=l,2,...cx)) 
occur in the final expression. This implies that the slowly varying (adiabatic) approximation 
Sj, where one merely replace E by E{t) in the formulas for pair production appears to be 
exact after Bose Enhancement effects are included. 
Now we present a derivation of eq. ([5]). 

The vacuum to vacuum transition amplitude in the presence of background field A is 
given by 

/|#t][rf^]e'/<i'-*'"I»l* ' ' ^ ' 

where for scalar field theory 

M[A] = {p- eAf - m\ with = (7) 

The one loop effective action is 

S^^^ = iTrln[(p - eAf - m^] - iTrln[p2 - m\ (8) 

where the Tr is given by 

TrO = j d^x< x\0\x > . (9) 



Since it is convenient to work with the trace of the exponential we replace the logarithm by 

Jo s 

Hence we get from eq. ([8]) 

Jo S 



We assume that the time dependent electric field E(t) is along the z-axis. We choose the 
Axial gauge As = so that only 

Ao = -E{t)z. (12) 
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is non-vanishing. Using eq. (fT2l) in eq. (fTTj) we obtain 

ds /'+°° r+oo r+oo r+oo 



t — / at < t\ dx < x\ dy < y\ dz < 

Jo S J —oo J —CO J —oo J —oo 



Inserting complete set of \pt > states / ci^prbr >< Pt\ = 1 we find (we use the normal- 
ization < q\p >= e*^'') 

^(1) = -li^ [ d^XT I d^pTe''^PT+-^"+'^)[r" dt < t\ dz < z\ 

{2'kY Jo s J J J-oo J~oo 

1 



is[{~i±+eE{t)zf~p^,] 



z>\t> - I dt I dz—]. (14) 



At this stage we use the shift theorem 



+00 



dy / dxfiiy) < x\ e-[("(^)"+'^^)'+''(^)+'=(^^)l \x > /sd/) 



oo 



= f dy rdxMy) <x-^M e-['^^(^)-^+^(i^H^(^)l \x - > Uy) 

J J-oo a[y) dy a{y) dy 

(15) 

to obtain: 

^(1) ' ds 



i2n) 



Jo S J J 



(16) 

where the z integration must be performed from — oo to +oo for the shift theorem to be 
applicable ^\. 

Inserting complete set of \pz > states (using / dpz\pz >< Pz\ = ^) as appropriate we find 

/r ]^ r+oo r r r+oo ^ ^ 

4^ - loo < 7 ^^^7 ^^^loo < ' + mdt^''' ^ 

< p.|e-[-^^^W^^+^^]|g, X q,\z + > \t >] 

• f (17) 



(27r)2 Jo s 



where 



F = /^°°^^ <t\J dPzJ dq.J^'^dz e^^P^e"^^^^ < p,|e^'[-"'^'W"'+^'] |g, > 
e^^^^e~*^^^|t > . (18) 
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It can be seen that the exponential e contains the derivative -n which operates on 



1 d „ 1 



< p^le*"*' ^ ^ ~^^''^\qz > "^'"^^ hence we can not move e ^(*) '''^^ to right. 

We insert complete sets of 1^; > states and \po >states using the completeness relations: 

/CXD /"OO 
dz\z><z\; 1= dpo\po><Po\ (19) 
-oo J —oo 



to find 



^ f27ry i J J ° i ° J ° i i i i i 



(27 



n+OO 



J dpo J dpQ j dpQ j dpQ j dpz j dq^ j dzi J dz2 J dz 



(27r)= 

^^tpo^izp^ ^ pQ|e-^^P-|p; > e-'^'P^ < p'o\ < zi\e''^-''^'^'^''+^'^\z2 > \po > 
e"^"^ < Pole^^'^'lPo > e-"«^e-**Po". (20) 

The matrix element, < |e*^[^^^^^(*)^^+^^l | Z2> can be evaluated in terms of the eigenstates 
of the Harmonic Oscillator with imaginary Q] and time dependent frequency uj{t) = 
ieE{t). The Hamiltonian 

H=pl + cu'{t)z' (21) 

can be diagonalized at every time t in terms of squeezed state time dependent creation and 
annihilation operators b^{t) and b{t) j^. Explicitly 

H = uj{t)[2b'^{t)b{t) + l]. (22) 

The eigenvalues of the Hamiltonian are just 

A„ = w(t)[2n + 1], n = 0, 1,2, (23) 

and the normalized eigenstates < z\nt > are the usual harmonic oscillator ones with uj — *■ 
uit). 

[pl + co\t)z'']\nt >= {2n + l)uj{t) \nt > (24) 
< z\n, >= Mz) = (^)^/'(^^^n(^V^)e-^^^ (25) 



Note that the above "instantaneous" wave functions with the eigen values [{2n + l)co'(t)] 
can also be directly obtained by constructing creation and anihillation operators from z, 
Pz and 00 (t) (like constant 00 problem) irrespective of the works in j^, [gI. By using the 
orthonormality relation of Hermite polynomials we find 

dz\M^)f = j dz\< nt\z >? = j dz(^f'^j^^\Hn{z^))\^e--^'^^" = 1 (26) 



and 



J dzijj^{z)ipm{z) = J dz < nt\z >< z\mt > 



rfz(^^)i/2^i/:(;.V^)//„(^V^)e--W^^ = S^m. (27) 
By using J dz\z >< z\ = 1 in the above equation we find 

<nt\mt>=6nm- (28) 

The completeness condition is 

Y,\nt><nt\ = l. (29) 

n 

Inserting complete set of harmonic oscillator states \nt > ( J2n l^t >< = 1) in eq. fl20|) 
we find 

F = XI XI / dt J dpo J dp'o J dpQ J dp^ J dp^ J dq^ J dzi J dz2 J dz 

— j dt J dpo J dpQ J dpo J dpo J dp^ J dq^ J dzi J dz2 j dz 



(2vr)= 



— J dt J dpo J dpQ J dpo J dpo J dp^ J dq^ J dzi J dz2 J dz 



(27r)3 „ J-00 

^itpo^iz{p,~g.) <p^|g-Etl)>.|p'^ > g-i^iP. < ^^1^^ > g-«ei?{t)(2n+l) ^ ^^|^^ > |^// > 

e'"^'?^ < p'o' I ^''^ > e-''Po' . (30) 

As advocated earlier the z integration must be performed from —00 to +00 for the shift 

theorem to be applicable ^. Performing the z integration we get 

r+00 . 

/ dze^'^P^-'^^^ = 27iS{pz - g.) (31) 
J —00 



and we find 

F = f^dt J dpo J dp', J dplj dp'^ I dp, I dz, I rfz2e^*^° <po|e-^^^1p'o > 

(32) 

Inserting further complete sets of \t > states we obtain 

F = j^-:^ Yl j dt J dti J dt2 J dpo J dp'f) J dp, J dp, J dp, J dzi J dz2e'^P° 
< j9o|e~^^^^|p'o > e'^'^P^ < p'olti >< til < zi\nt > e-^^^W(2n+i) ^ ^^|^^ > 
\t2 >< t2\Po > e'^'"^ < Pole^^^^p'o" > e~**< 

= H j dt J dti J dt2 J dpo J dpQ J dpi J dp^ J dp, J dzi J dz2 

1^2 > e^*2Poe^^2P. < p'^\e^ip^ \p'^' > g-^H". (33) 

Since < nt\z > does not contain any time derivative ^ operator (see eq. (125|1 ) we find (by 
using < ti\t2 >= S{ti - t2)) 

F = ^ y dt J dti J dpo J dpQ J dpi J dp^ J dp, J dzi J dz2 

e'^'P^ < p'ole^^'^'lPo > e-'*<. (34) 

1 d 1 d 

Since < po\e * bo > and < Po|e^('' \Po > are independent of the time derivative op- 
erator ^ (they only depend on c-numbers po, p'q, p'q and p'^) we can take < PqIc"^^^"^ |pq" > 
and e~**^o" to the left. We find 

F = -j^-;^ ^ j dt j dti j dpo J dpi J dpi J dpi' J dp, J dzi J dz2 

^it{po~p'i') < pl\e^)iiP^\p';; >< p^\e'wr)iiP^\p'^ > e--iP^e-**^^'o < zi\nt, > e-^^^(*i)(2n+i) 
<nt,\z2> e'''Poe'''P\ (35) 

1 d 1 d 

It can also be shown that < po|e ''^ |Po > aiid < Po|e^(*' |po' > are independent of t. 
This can be shown as follows: 

< pol/(t)|b'o >= / dt'J dt" j dp';;' < p,\t' X t'lmit" X t'>- X p'i;'\j^\pi > 



= J dt' J dt" J dp'^'e-'''P' 5{t' - t")f{t'y>o" ip'^ 5{p'^'-p'o) = ip'o j dt' e-**'(^'°-^'o)/(t') 

(36) 

which is independent of t. Hence we can easily integrate over t in eq. ( l35l) . 

Integrating over t in eq. fl35|) by using dte^*^^^°~^o'^ = 2n6{po — p^), we find 

^ = / dti J dpo J dp'oj dplj dp, J dz^ j dz2 < Pole^^P^o ><Po|e"^^^1Po > 

Since / dpo\po >< po\ = 1 and / rfp^e*^^*^^^"^^-' = 2n6{zi — Z2) we find 

F = / dtij dp'^j dp^J dz,J dz^e-'^^P^ < z,\nt, > e— < nt,\z2 > e 

= J-y^ [ dt [ dp', [ dz\ < z\nt > p e-^^^W(2n+i)_ 
(2vr) n J J J 

Since the harmonic oscillator wave function is normalized (see eq. (l26l) ): 

j dz\ < z\nt > P = 1 (39) 

we find 

^ - E / * / *oe-«-<'..-.) ^ / * / *05-i^ (40) 

Using the Lorentz force equation 

dp^ = eF^^dx" (41) 
we find (when E{t) is along z-axis, eq. ffT^ ) 

dpo = eE{t)dz. (42) 

Hence 

dt f dz-^^^—. (43) 



(47r) 7 7 sinh(seE(t)) 
Using the above expression for F in eq. (ITTl) the effective action becomes 

^(1) 



fd^:, [ ^2 ]. (44) 

167r3io s J J ^ smh{seE{t))^ ^ ' 
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The series expansion for is ffiven by 



+ 2^E4A:^- (45) 



sinhx X ~^ TT^n^ + x'^ 
Hence 

The s-contour integration is straight forward . Performing the s-contour integration 

around the pole s = j^^pyy we find for the probabihty W of pair production (which is twice 
the imaginary part of the effective action) 



W = 2Im5(^) = ^ E ^ / / rfV|e^(t)|e-""W^. (47) 

8^ n=i n J J 



Therefore the probabihty of producing a charged particle per unit volume per unit time with 
transverse momentum pt is given by 

- ' ^ -Log[l + e i-^wi ] (48) 



d'^xd?pT Stt^ 

which reproduces eq. ([5]). 

To conclude we have studied, for the first time, the Schwinger mechanism in the presence 
of arbitrary time dependent background electric field E(t) by directly evaluating the path 
integral. We have obtained an exact one-loop non-perturbative result for the probability of 
charged scalar particle-antiparticle pair production per unit time per unit volume per unit 
transverse momentum (of the particle or antiparticle) from the arbitrary time dependent 
electric field E{t). We have found that the exact non-perturbative result is independent of 
all the time derivatives where n = l,2,....oo. 

The quantum Back Reaction problem for the Electric Field is studied in sl- For the Back- 
Reaction problem, the equation for the charged scalar field is supplemented by the Back- 
Reaction equation for the time evolution of the Electric Field , which is the semiclassical 
Maxwell's equation: 

^ -(j'W) (49) 



dt 

where j^(t) is the induced (quantum) current coming from the produced charged scalar 
mesons, and the expectation value is taken with respect to the initial adiabatic vacuum 
state of the scalar field theory. The exact numerical solution of the coupled semiclassical 



Maxwell equation and the equation for the time evolution of the quantum field were 
compared in those papers to the solution of of the semi-classical transport equation having 
a source term for pair production which utilizes Schwinger's constant electric field E result 
but with the replacement E —>■ E{t). The transport approach agreed with a coarse graining 
of the quantum field theory solution. This result seems to be consistent with our finding here. 
However in those simulations the Schwinger source term was modified by Bose Enhancement 
effects for boson pair production and Pauli-blocking effects for fermion pair production 
js, [o]. These enhancements or (Pauli-blocking) effects are not present in our calculation 
and need to be better understood. Particle production from arbitrary time (and space) 
dependent background fields is important in early universe dynamics and in the production 
and equilibration of the quark-gluon plasma at RHIC and LHC 
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